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ABSTRACT 

The distribution of eigenvalues of N x N random matrices in the limit N — > oo is the 
solution to a variational principle that determines the ground state energy of a confined 
fluid of classical unit charges. This fact is a consequence of a more general theorem, 
proven here, in the statistical mechanics of unstable interactions. Our result establishes 
the eigenvalue density of some ensembles of random matrices which were not covered by 
previous theorems. 
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I. INTRODUCTION 

Since the pioneering work of Wigner [37,38,39], there has been a considerable effort to 
understand the statistics of eigenvalues of N x N random matrices. The problem has three 
scales: (i) the density of eigenvalues which converges to a deterministic limit as N — > oo; 

(ii) the fluctuations of order one around this deterministic density; (iii) the local statistics 
on the scale of the typical distance between eigenvalues. Item (i) depends on the particular 
matrix ensemble while (ii) and (iii) are "universal," in the sense that they depend only 
on some overall matrix characteristics (e.g. the matrices being real and symmetric). The 
"classical" results are reviewed in [28]; see also [31] for a collection of early work. Recently, 
in the context of the double scaling limit of 2D quantum gravity [10], (ii) and (iii) have 
been studied at the edge of the support of the density of states where novel universality 
classes occur. Among recent work on (ii) we also mention [29,34], and [30,4,9] regarding 

(iii) . In our paper we will consider only the largest scale (i). 

For various N x N random matrix ensembles of the form 

m^ N \dM) = Q(iV)- 1 e- KiVTt ^ M ^M, (1.1) 

the joint probability distribution for the N eigenvalues Ai, ...,\n (which may be real or 
complex) is identical to the (configurational) canonical ensemble at inverse temperature 
P of N unit point charges at positions Ai, Ajv G A C R 2 . The region A can be all of 
1R 2 , the unit disk Bi C M 2 , the unit circle S 1 , the entire real line R, or some other set, 
depending on the type of random matrices. This joint probability distribution has the 
general form 

=Q( N \f3)- 1 exp(-f3H ( > N) )d\ 1 ---d\ N (1.2) 

on A N , where d\k is the uniform measure on A and Q^ N \j3) the normalizing partition 
function. The classical Hamiltonian, H^ N \ is of the form 

H^ N \X U ...,X N ) = G(\j,\ k ) + J2 F(^k)+NV(X k ) 7 (1.3) 

l<j<k<N l<k<N 

where G(Xj, A&) = G(Afc, Xj) is (2ttx) a Green's function for —A in two dimensions, F(X) = 
lim ?? ^A(G r (A, rj) + In |A — r)\) is the regular part of G, and V(X) given in (1.1). 
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Let us list a few examples. The joint eigenvalue distribution, fi^ N \ of the Gaussian 
ensembles (i.e., V(M) = M^M in (1.1)) for the real symmetric [22], the complex Hermitian 
[14] , the general complex [20] , and the Hermitian self-dual quaternionic [14] N xN random 
matrices takes the form 

d " W = 0<4sT I! Hi-M' II e-^'dX,; (1.4) 

^ l<i<j<N l<k<N 

see also [28]. Clearly, this corresponds to (1.2), (1.3), where G(Ai,A 2 ) = — ln|Ai — A 2 | is 
the free space Green's function (whence F = 0) and V(X) = |A| 2 a quadratic potential. The 
eigenvalues for real symmetric, complex Hermitian, and Hermitian self-dual quaternionic 
matrices are real. Therefore, the charges are confined to R, i.e. A = R. The parameters 
have the values ,3 = 1,2,4 and k = 1,2,2, respectively. For each of the associated ensembles 
of unitary matrices the charges are confined to the unit circle, i.e. A = S 1 . Then f3 is 
unmodified but k = [13]. For general complex random matrices, the eigenvalues are 
complex, corresponding to unconfined charges, i.e. A = R 2 = C, and (3 = 2, k = 2 [20]. 
The joint eigenvalue distribution for the general real quaternionic N x N matrices takes 
the more complicated form [20] 

^ iN) = oiWF) II lA.-A/lA.-A*^ J] \X k -Kfe~^ 2 dX k (1.5) 

with /3 = 2,k = 1,A = IR 2 = C. Although apparently not noticed previously, (1.5) 
can be interpreted as a configurational canonical Coulomb ensemble, with a Hamiltonian 
(1.3) in which now G(Ai, A 2 ) = — In |Ai — A 2 | — In |Ai — | is the Green's function of —A 
for the upper half space R x R + = C + equipped with a perfectly dia-electric condition 
at its boundary dC + (the real axis) and asymptotic free conditions at infinity, extended 
symmetrically to R 2 , having a regular part given by F(X) = — ln(2|Jm(A)|). Moreover, 
V(X) = |A| 2 . This Hamiltonian describes Coulomb point charges interacting via the free 
space Green's function — ln|A — A'| amongst each other and also with A^ identical image 
charges with respect to the line Xm(A) = 0. Since the interaction of a charge with its 
own image contributes only an amount F/2 to the Hamiltonian, F/2 + NV is now to be 
counted as the external potential. 
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For symmetric random matrices, properties of the eigenvalue statistics have been 
computed in great detail for arbitrary N, using explicit expansion techniques [28, 27, 20], 
group theoretical methods [13], the method of orthogonal polynomials [28,30,4,9], as well 
as some more recent developments in soliton theory and two-dimensional quantum gravity 
[1,2]. The beautiful connection to two-dimensional Coulomb systems suggests to use the 
general methods of statistical mechanics when the above algebraic methods fail. Even in 
the exactly solvable situations the statistical mechanics approach may provide us rather 
readily with certain relevant asymptotic (N — > oo) results which to extract from the exact 
finite N solutions would require quite tedious and lengthy computations. Furthermore, as 
Dyson has pointed out [13], in the framework of statistical mechanics the limit N — > oo 
makes sense for arbitrary (3 and not only for the discrete values of j3 = 1, 2, 4. Thereby one 
achieves a "thermodynamic view point" which yields valuable new insights into random 
matrices. Early results in this direction are in [13,40], and more recent ones in [8,17,18]. 

The prime example of the exploitation of the Coulomb analogy is Wigner's [39] elec- 
trostatic derivation of his semi-circle law 



for the eigenvalue density p(X) in the limit N = oo of the Gaussian ensemble (1.1) in 
the case of real symmetric random matrices. Wigner [39] argued, heuristically, that when 
N — > oo, the eigenvalue density for (1.4) can be obtained from a variational principle for 
a continuum charge density p(X) of total charge 1, restricted to the real line, that satisfies 
the requirement of mechanical force balance between its own electrostatic force field and 
the applied force field — <9a|A| 2 = — 2A. Previously [37] he had proved, by the method of 
moments, that (1.6) holds true for a Bernoulli ensemble of bordered random sign matrices, 
which suggested that (1.6) was the limiting law under more general circumstances [39]. In 
[38,39] he announced that (1.6) can indeed be proved to hold for a wider class of ensembles 
under a mild set of conditions, including the Gaussian real symmetric ensembles. 

Interestingly in itself, entropy plays no role in Wigner's variational principle, which 
is concerned only with the ground state energy of the classical continuum Coulomb fluid. 
To get an intuitive idea how this can arise from the canonical measure (1.4) with fixed (3 




(1.6) 
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(< oo), we rewrite (1.4) as 



AT 

l<i<j<N k=l 



n dx k- (i-7) 



Kk<N 



Apparently, the limit A" — > oo for (1.7) is a simultaneous thermodynamic anc? zero- 
temperature limit for an unstable Hamiltonian with mean-field scaling. If in (1.7) we 
replace (3N by (3N , with Ao fixed, and then let A" — > oo, we obtain the variational prin- 
ciple of [6,23] for a continuum free energy functional at inverse mean-field temperature 
Pmf = /3Ao, see also [26]. Letting Ao — > oo subsequently, the entropy contribution to the 
free energy drops out, giving formally Wigner's variational principle. With a leap of faith 
one may thus expect that the limit AT — > oo in (1.7) will give the same result directly. 

Recently, Boutet de Monvel, Pastur, and Shcherbina [5] have studied the limit A" — > oo 
of the joint eigenvalue distribution of real symmetric random matrix ensembles (1.1) for 
a large class of V(M), satisfying certain regularity conditions. They prove that in the 
limit A" — > oo the n-th marginal measure of (1.4), with A 2 replaced by V(A), factors 
into an n-fold tensor product of identical one-particle measures whose density is precisely 
that of the two-dimensional Coulomb fluid restricted to a line, satisfying the equations of 
electrostatic equilibrium in the applied potential field V. The proof in [5] is based on the 
classical Stieltjes transform. It is carried out explicitly for (3 = 1, but the method covers all 
(3 > (in particular, including (3 = 2,4), as well as certain Holder continuous many-body 
interactions. 

We here generalize this result in [5] to a wider class of interactions and an arbitrary 
(finite) number of space dimension. Our result covers also the limit A" — > oo of the 
eigenvalue distribution of Ginibre's complex and real quaternionic Gaussian random matrix 
ensembles [20]. However, our method is very different from that in [5]. Instead of using 
the Stieltjes transform, we adapt the strategy of Messer and Spohn [26], Kiessling [23] 
and Caglioti et al.[6] to the combined mean-field and zero temperature limit. Interestingly 
enough, we can work without the detailed control of [5] on the finite- A" marginal densities, 
and in this sense our proof of the variational principle is also considerably shorter and 
simpler than the proof in [5]. 
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II. MAIN RESULT 

We now prepare the statements of our main result. Let A C t d be closed and con- 
nected, and let dx denote uniform measure on A. Note that A may be all of M d , or it may 
be a lower dimensional manifold, e.g. the sphere S d_1 . We denote by P(A) the probability 
measures on A, and by P S (A N ) the permutation symmetric probability measures on the 
infinite Cartesian product A N . We recall the decomposition theorem of de Finetti [16] and 
Dynkin [12] (see also [21,15]), which states that \i G P S (A N ) is uniquely presentable as a 
linear convex superposition of product measures, i.e., for each fx G _P S (A N ) there exists a 
unique probability measure v(dg\fi) on P(A), such that for each nGN, 

fi n (d n x)= [ v(d g \,u)Q® n (d n x), (2.1) 

JP(A) 

where fi n denotes the n-th marginal measure of fi. This is also the extremal decomposition 
for the convex set P S (A N ), see [21]. 

To establish the limit N — > oo for measures of the form (1.2), (1.3) only some general 
properties of the interactions enter. Also, the particular value of (3 plays no role, and 
it is convenient to absorb it into the Hamiltonian. Therefore, we study the sequence of 
probability measures 

^ JV )(^) = -^ y exp(-^)(x 1 ,..,x iV )) J] dx k (2.2) 

^ l<k<N 

for Hamiltonians of the form 

#W( Xl ,..., Xi v)= Yl «^( x ^ x j)+ Yl u(^)+Nv(^). (2.3) 

l<i<j<7V l<i<N 

The pair interaction w and one-particle potentials u, v satisfy the following conditions. 
Condition on ty(x, y): 

(CI) Symmetry : w(x, y) = w(y,x) 

Condition on U(x) = -u(x) + w(x): 

(C2) Integrability : e~ u{ ^ G L X (A, dx) 
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Conditions on W(x, y) = w(x, y) + v (x) + v(y): 

(C3) Lower semicontinuity : W is l.s.c. on A x A 



(C4) Integrability : W(x, y) G L 1 (A 2 , e"^ x )da; <g> e"^ y )dy) 

In case of an unbounded A we also need a growth condition at infinity. Let VF_(x) = 
min y W(x, y). 

(C5) Confinement : lim VF_(x) = oo , uniformly in x 

THEOREM: Let A C M d be closed and connected, and letw, v and u satisfy the conditions 
(CI) — (C5). Consider (2.2) as extended to a probability on A N . Then there exists a 
[i G P S (A N ) such that, after extraction of a subsequence ^ N \ 

lim = ^ . (2.4) 

For each limit point \x, the decomposition measure u(dg\n) is concentrated on the subset of 
-P(A) which consists of the probability measures g that minimize the functional 

S(g) = \g®\W) (2.5) 

over -P(A) . 

In general we have little information on the decomposition measure v(dg\/i). More is 
known for regular mean-field Hamiltonians, see [25]. However, if it can be shown, as is the 
case for many random matrix ensembles, that (2.5) has a unique minimizer, say go, then 
in (2.4) we have in fact convergence, and the limit is of the form 

\x n = QT- (2-6) 

As discussed in [35], the factorization property (2.6) is equivalent to a law of large 
numbers. Consider the eigenvalues averaged over some continuous test function /, 

N 

s 



1 N 



Then (2.6) implies that, for all such /, 



lim 



(f) N = [ Qo(d\)f(X) 
J A 



(2.8) 



in probability. We summarize as 

COROLLARY: // (2.5) has a unique minimizer, go, then the weak law of large numbers 
(2.8) holds for all f G C°(A). 

Applications of our theorem to random matrix ensembles are presented in the con- 
cluding section. 

III. PROOF OF THE THEOREM 

The proof looks somewhat technical because we work under fairly minimal assump- 
tions on W. But, in essence, we only have to show, through sharp upper and lower bounds, 
that the pair-specific free energy converges to the minimum continuum energy. The re- 
mainder of the theorem follows from the permutation invariance of the measures. 

We define the absolutely continuous (w.r.t. dx) a-priori probability measure on A, 



Hoidx) = Z^e-u^dx. 



(3.1) 



For each g( N ) G P(A N ) its entropy w.r.t. [i® N is defined by 




(3.2) 



if q( n > is absolutely continuous w.r.t. a-priori measure nf N , and provided the integral in 
(3.2) exists. In all other cases, (g {N) ) = -oo. 



LEMMA 1: The relative entropy (3.2) is non-positive, 



SW(gW) <0. 



(3.3) 



Proof of Lemma 1: 



A standard convexity argument, see [15]. 



□ 
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We introduce the symmetric Hamiltonian 

^(xi,.,x,) = i W ^j^)- ( 3 - 4 ) 

^<j^k<N 

For qW e P(A N ), (px) its Helmholtz free energy (or just free energy) is now defined by 

jrW = e W (#)) _ 5 W (3 _ 5) 

if the right side exists, and by (g( N >) = +00 in all other cases. 

LEMMA 2: The free energy (3.5) takes its unique minimum at the probability measure 

^ N) {d N x) = ^ y exp(-KW(x 1 ,...,x Ar ))^(A), (3.6) 

where = QW/Ztf. Thus, 

min J* N )M N ))=FWU N )) = -hiZW. (3.7) 

Proof of Lemma 2: 

The variational principle is verified by a standard convexity estimate [15,32], which 

shows that J=W (qW) -J 7 ^ (fj,W) > 0, with equality holding if and only if qW = 

The second identity in (3.7) is verified by explicit calculation. □ 

Notice that the canonical probability measure (3.6) is just (2.2) rewritten in terms of 
an d ^ 

LEMMA 3: The functional £(g) has a finite minimum over -P(A). 

Proof of Lemma 3: 

For bounded A the lower semicontinuity (C3) establishes that £(g) is a lower semi- 
continuous functional for the topology of measures in -P(A). For an unbounded region, the 
same conclusion holds because of (C3) and (C5). Now the claim follows from standard 
facts about lower semicontinuous functionals [32]. □ 



In the following, let go denote a minimizing probability measure for (2.5), and let 
E = £(g ). 
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LEMMA 4: The pair specific free energy is bounded above by 



lim sup 



< E . 



(3.8) 



We will prove Lemma 4 in two blocks, distinguishing minimizers go with finite entropy 
from those with negative infinite entropy. For instance, Wigner's semicircle density (1.6) 
has finite entropy, but other ensembles can have a singular measure as eigenvalue "density." 

Proof of Lemma 4 

Assume first that the entropy of go is finite, i.e. (recalling Lemma 1), 



Notice that the non-positive constant So is ^-independent. We then can use Lemma 2 to 
estimate for all iV that 

-N~ 2 \nZ^ <N~ 2 ^ N \gf N ) = (l-N- 1 )£(g )-N- 1 S^(g ) . (3.10) 

The identity in (3.10) follows from (CI). By (3.9), (qq) = Sq, and So is finite and 
iV-independent, and by Lemma 3, £{qq) = E is finite and iV-independent. Thus (3.8) in 
the finite entropy case follows by taking N — > oo in (3.10). 

Assume now that (3.9) is false, so that S^(go) = — oo. In that case, (3.10) becomes 
useless. Now, since the functions are dense in -P(A), the obvious way out is to modify 
the above argument and to work with a regular approximation to go- However, since S(g) 
is only lower semicontinuous, we have to employ also a continuous approximation to S(g). 

By (C3) for bounded A, and by (C3), (C5) in case of an unbounded A, W is the 
pointwise upper limit of a continuous increasing map 7 1— > W 1 G C°(A x A) that is 
uniformly bounded below, see [32]. By (CI), we can assume that W 7 (x, y) = W 7 (y, x). In 
the following, let Kj be defined by (3.4) with W replaced by W 7 , and let fij N ^ (d N x) and 
Zj N \ respectively (g^), be defined by (3.6), respectively (3.5), with replaced 



(qo) = S e( 



oo,0] . 



(3.9) 
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Since W satisfies (C3), and W 1 is of class C°, and W 1 / W pointwise, for each 
positive e<C 1 we can find a r y e such that, simultaneously, 

Eo ~ \qT^) < \e (3.11) 

and 

limsup (-N- 2 \nZ {N A < limsup (-N~ 2 In z!f A + (3.12) 

whenever 7 > 7 e . Moreover, since W 7 is bounded below on A x A uniformly in 7, and 
W 7 G C°(A x A), for each 7 e we can find a measure g e G -P(A) that is equivalent to a 
positive function of class Cq°(A), such that 



^? 2 (w 7 .)-i^ 2 (w 7 .) 



<^e. (3.13) 



On the other hand, given any 7 and any g s G -P(A) of class C^°, we have the estimate 

- In Zf = ^ (#)) = min ) 

< {gf N )=N(N-l) 1 -gf 2 (W 1 )+NS^(g s ), (3.14) 

where the first line is the analog of Lemma 2, the inequality obvious, and the last identity 
an explicit computation, using the symmetry of the W 1 . In particular, for any given e 
we can choose 7 = 7 e and g s = g e in (3.14), then multiply (3.14) by A^ -2 and take the 
limsup. Since g e is of class Cq°, we have \S^(g e )\ = C(e) < 00, independent of N, 
whence iV -1 S (1) (e e ) -> 0. Next we use (3.11), (3.12), (3.13), and the triangle inequality, 
and conclude 

limsup (-AT- 2 In Z {N) ) < \ e f 2 {W lt ) + \e < E + e , (3.15) 

N^oo ^ / 2 3 

for arbitrarily small e. This proves (3.8) for the infinite entropy case. 

The proof of Lemma 4 is complete. □ 

LEMMA 5: The sequence N 1— > /j( N ^ given by (3.6) is compact for bounded A, and tight 
for unbounded A. 
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Proof of Lemma 5: 

If A is bounded then it is also compact, for A C IR d is closed, and in that case -P(A) 
is compact for the topology of measures. By Tychonov's theorem, P S (A N ) is now compact 
in the product topology. Hence, for bounded A, the sequence (3.6) is compact. 

If A is unbounded, we have to estimate the contribution from outside of to the mass 
of the n-th marginal fi n N ^ of (3.6), when N oo. Recall that a sequence of probability 
measures fi n N ^ is tight if for each eC 1 there exists a R = R(e) such that (Br) > 1 — e, 
independent of N, see [11]. Since our marginal measures are compatible (i.e., fi n (d n x) = 
/J,m\d n x £g> A m_n ) for m > n) and permutation symmetric by (CI), it suffices to prove 
tightness for any particular n. We pick n = 2. 

We notice that by (C3) and (C5), 

min W(x, y) = W > -oo . (3.16) 

(x,y)£AxA 

Since (3.6) is invariant under the transformation W — > W + C, we can even assume, 
without loss of generality, that W > 0. We then have the following sandwich bounds, 
independent of N, 

0<^ 2 N \w)<^ 2 (W), (3.17) 

with n® 2 (W) < oo, by (C4). The lower bound in (3.17) is obvious, for W > 0. To prove 
the upper bound in (3.17) we use the strategy of [23]. We can replace by otK^ 

in (3.6), with < a < oo, so that -2A^" 2 (1 - jV _1 ) lnZ^ = T N (a) is now a function 
of a. Clearly, rAr(O) = 0, and W > implies rAr(a) > as well as T' N (a) > 0, while 
the Cauchy-Schwarz inequality implies F'^(a) < 0. Moreover, Jensen's inequality, applied 
w.r.t. (j,q N , and (C4) imply 

T N (a) < afif 2 (W) . (3.18) 

Obviously (j,q 2 (W) is A^-independent. Thus, Tat (a) is a nonnegative, increasing, concave 
real function, bounded above by (3.18), and mapping zero into itself. A simple geometrical 
argument now reveals that the slope of any tangent to rV(a:) never exceeds the slope of 
the ray on the r.h.s. of (3.18), i.e., T' N (a) < ^ 2 {W). But 1^(1) = (J& N) (W), which proves 
the right inequality in (3.17). 
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Now pick eCl arbitrary. By (C5) we can find a R = R(e) such that 

inf %y)>V(f). (3.19) 

(x,y) B| C 

Let x denote the characteristic function of the complement of B# in A. We then have the 
chain of estimates 

> inf W(^y)^\ X ^) > \p$\W) (l-»i N) (B£)) . (3.20) 

(x,y) g B 2 R e V / 

Division of (3.20) by e _1 /j,f 2 (W) and a simple rewriting reveals that, independent of N, 

»l N) (B£)>l-e, (3.21) 
which was to be shown. The proof is complete. □ 
LEMMA 6: The pair specific free energy is bounded below by 

lrniinf lnZ (Ar) ) > E . (3.22) 

Proof of Lemma 6: 

By Lemma 1, (^) < 0. Therefore, 

-]nZW>nW(KW). (3.23) 

By (CI), 

^ {N K^ N) )={l-^) 1 r i N) (W). (3.24) 

Now pick a converging subsequence of (3.6), P S (A N ). Such a converging 

subsequence exists by Lemma 5 and the Bolzano- Weierstrass theorem. Then, by (C3), we 
have 

liminf^f \W) >»2(W), (3.25) 
while 1 - N'' 1 -> 1 trivially. Thus, 

lrniinf (-N" 2 ln(Z^)) > -^(W) . (3.26) 
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Finally, using the representation (2.1), we see that 



\»2{W)= [ u(d 6 \fj)£(Q)>£(Qo), 



(3.27) 



and the proof of Lemma 6 is complete. 



□ 



Proof of the Theorem. 



By Lemma 4 and Lemma 6, 



lim {-N~ 2 \nZ {N) ) = E . 



(3.28) 



Recalling (3.26) and (3.27), we see that (3.28) implies 




for every limit point fi of fj,( N K Equation (3.29) in turn implies that the decomposition 



measure u(dg\fx) is concentrated on the minimizers of £(g), for assume not, then 



We are now also in the position to vindicate our remark on the existence of the limit 
in (2.4) in case the minimizer go is unique. Indeed, in that case the set of limit points of 



IV. APPLICATIONS 

With the specifications in (2.2), (2.3) that x = A G A C M 2 , v(X) = /3V(X), u(X) = 
f3F(X), and w(X,rj) = (3G{\ 77), where G is a Green's function for —A in 2D and F 
its regular part, our theorem characterizes the limit N = 00 of (1.2), (1.3), which for 
(3 = 1, 2,4 is the joint eigenvalue distribution of various random matrix ensembles of the 
form (1.1). The decomposition measure of the limit is concentrated on the ground state(s) 
of the electrostatic energy functional e(g) = (3~ 1 £(g) of a charged continuum fluid with 
"charge density" dg/dX (which may be a singular measure) of total charge 1, subject to 
an external potential V. Explicitly, the energy functional reads 




which contradicts (3.29). The proof of the Theorem is complete. 



□ 



{^ N \ N = 1, 2, . . .} consists of the single measure. 



e(g) = -Q® 2 (G) + g(V). 



(4.1) 
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The regular part, F, of G does not contribute to the limit. We list a few examples. 

IVa. Real symmetric, complex Hermitian, and quaternionic self-dual Hermitian matrices 
As mentioned in the introduction we have A = R, G(X, rj) = — In j A — 77 j , (3 = 1,2, 4, 
and k = 1,2,2, respectively. Our electrostatic variational principle (VP) for (4.1) then 
becomes the VP of Boutet de Monvel et al. [5], but here with a slightly wider class of 
potentials V. In particular, for (3 = 1 we can allow continuous V with V(X) ~ (1 + e) In |A| 
asymptotically, as compared to Holder continuous V with V(X) ~ (2 + e) In |A| in [5]. The 
VP has been studied extensively in [33]. A unique minimizer is known to exist under 
certain regularity conditions on V. Of course, for V(X) = |A| 2 , the quadratic potential of 
the Gaussian ensembles, the minimizer of (4.1) is given by Wigner's semicircle law (1.6). 

IVb. General complex matrices 

We have A = R 2 , G(X, 77) = — In j A — 77 1 . In this case our variational principle for (4.1) 
generalizes the VP of [5] to two-dimensional domains. Under mild conditions on V, and 
in particular for all our examples, it can be shown that the minimizer is unique. 

We consider only the Gaussian ensemble with kV(M) = M^M in (1.1), whence 
V(X) = I A 1 2 / 2 in (1.3), and (3 = 2 in (1.2). The minimizer of (4.1) is given by 

dg = 7t- 1 XBi (X)dX, (4.2) 

where % B (A) is the indicator function of the unit disk Bi in R 2 . This result can also be 
obtained from Ginibre's exact finite N formula, see [20]. 

IVc. Complex normal matrices 

We have A = R 2 , G(X,rj) = -ln|A - rj\, and (3 = 2. Consider first (1.1) with 
kV(M) = ln(l + M^M) 1+1 / N . Then in (1.3) we have V(A) = - ln^dAI)] 1 / 2 , where 
Pc(€) = 7r_1 (l + £ 2 ) -1 is the density of the Cauchy distribution, and F is replaced by V. 
With these identifications (C5) is violated, but (1.2) is well defined for all (3 > 1, and the 
minimizer of the electrostatic energy functional is found to be 

^ = 7T" 1 (l + |A| 2 )" 2 rfA. (4.3) 

The measure (4.3) has geometrical significance. Recall that |J| 2 (A) = 4/(1 + |A| 2 ) 2 
is the Jacobian of the stereographic projection map S 2 — > R 2 , arranged such that the 
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equator of § 2 coincides with the unit circle in IR 2 . Therefore, (4.3) is the stereographic 
projection onto the Euclidean plane of the uniform probability measure on § 2 . Also the 
finite- N measure (1.2), with (1.3) specified as above, is itself a stereographic projection 
onto Euclidean space of a canonical ensemble measure of N point charges in the two-sphere 
S 2 . For (3 = 2, this spherical ensemble is given by (1.2) with Hamiltonian 

hW(\ 1 ,...,\ n ) = - HAj-Aifcl, ( 4 - 4 ) 

l<j<k<N 

where Xj G S 2 , \Xj — A&| is the chordal distance on S 2 , and dXk in (1.2) now means uniform 
measure on § 2 . (If f3 ^ 2, a non-constant one-particle potential has to be added to H^ N >). 
Our theorem applies directly to this ensemble on (S 2 ) x7V . We arrive at (4.3) by taking the 
limit N — > oo on the sphere, arguing that the minimizer is a constant, and projecting the 
result onto the Euclidean plane. 

In a sense, the spherical ensemble is the counterpart on S> 2 to Dyson's circular ensem- 
bles on S 1 , although these ensemble are related to random matrices in a different manner. 
The statistical mechanics of the spherical ensemble has been studied in some detail in [7,19], 
using exact algebraic techniques. The spherical ensemble is also related to J.J. Thomson's 
celebrated problem of determining the minimum energy configuration of iV point charges 
on S 2 [36], which recently has received much attention in physics [3], topology and Knot 
theory [24]. 

As a second example, consider entries restricted by HMH^ < 1 for every M, with 
uniform distribution otherwise. This corresponds to V(M) = lim t _ >00 (l + tanh[£(7 — 
MM*)]) -1 . Alternatively, we can choose A = Bi and V = in (1.2). The unique 
minimizer of the corresponding variational principle is the electrostatic charge distribution 
on a circular perfect conductor. It is well known that any surplus charge accumulates at 
the "surface," i.e., the minimizer is given by a Dirac mass concentrated uniformly on the 
boundary of the unit disk, 

dg = 5 § i(dx). (4.5) 
With the help of our corollary this gives us the following. 

Proposition: Let (f) N be defined as in (2.7), the summation running over the eigenval- 
ues, restricted to B\ cC, of a complex normal N x N random matrix whose free entries 
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are uniformly distributed otherwise. Then, in probability, (f) N — ► (27r) 1 J ^ / (e 1 ^) d</? 
as A" — > oo. 

This result is worth rephrasing in less technical terms. As far as averages over the 
spectrum are concerned, an infinite normal random matrix with eigenvalues in the unit 
disk and independent entries uniformly distributed otherwise, is almost surely equivalent 
to some infinite unitary matrix. 

IVd. Real quaternion matrices 

We consider only the Gaussian ensemble. The joint probability density is given by 
(1.5) with A = M 2 , j3 = 2, k = 1, [20,28]. The limiting eigenvalue distribution now 
minimizes the slightly more complicated electrostatic energy functional 

e(g) = \ e ® 2 (- In | A - v \ - In |A - v *\ + |A| 2 + M 2 ) , (4.6) 

where r]* is the mirror image of r\ with respect to the real axis. The minimizer of (4.6) is, 
once again, unique and given by the measure (4.2). 
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